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Abstract 

We develop a new unified theoretical approach enabling us to non-perturbatively study the effect of electron-electron 
interactions on weak localization in arbitrary arrays of quantum dots. Our model embraces (i) weakly disordered conductors 
(ii) strongly disordered conductors and (iii) metallic quantum dots. In all these cases at T — > the electron decoherence 
time is determined by the universal formula t v o ~ grr>/ \n(Ec/S), where g, rr}, Eq and S are respectively dimensionless 
conductance, dwell time, charging energy and level spacing of a single dot. In the case (i) this formula yields t v q oc D 3 / In D (D 
is the diffusion coefficient) and matches with our previous quasiclassical results [D.S. Golubev, A.D. Zaikin, Phys. Rev. Lett. 
81 (1998) 1074], while in the cases (ii) and (iii) it illustrates new physics not explored earlier. A detailed comparison between 
our theory and numerous experiments provides an overwhelming evidence that zero temperature electron decoherence in 
disordered conductors is universally caused by electron-electron interactions rather than by magnetic impurities. 
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1. Introduction 

Quantum interference of electrons in mesoscopic 
conductors manifests itself in a number of fundamen- 
tally important phenomena which can be directly 
observed in modern experiments. One of them is the 
phenomenon of weak localization (WL) [1,2,3]. In 
the absence of interactions electron wave functions 
preserve their coherence and, hence, quantum inter- 
ference remains efficient throughout a large part of the 
sample making WL a pronounced effect. Interactions 
between electrons and with other degrees of freedom 
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may limit phase coherence thereby making quantum 
interference of electrons possible only within a finite 
length scale L v . This so-called electron decoherence 
length as well as directly related to it decoherence 
time t v = L^,/D (where D is diffusion coefficient) are 
crucial parameters indicating importance of quantum 
effects in the system under consideration. 

At sufficiently high temperatures quantum behavior 
of electrons in disordered conductors is usually sup- 
pressed due to various types of interactions. However, 
as temperature gets lower, certain interaction mech- 
anisms either "freeze out" or become less efficient in 
destroying quantum coherence. As a result, both L v 
and t v usually grow with decreasing temperature and 
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quantum effects become progressively more important. 

Should one expect L v and t v to diverge in the 
limit T — > 0? While some authors tend to give a 
positive answer to this question, numerous experi- 
ments performed on virtually all kinds of disordered 
conductors and in all dimensions demonstrate just 
the opposite, i.e. that at low enough T both deco- 
herence length and time saturate to a constant and 
do not anymore grow if temperature decreases fur- 
ther. The list of corresponding structures and exper- 
iments, by far incomplete, includes quasi-ld metallic 
wires [4,5,6,7,8,9,10,11,12,13,14,15,16], quasi-ld semi- 
conductors [17,18,19], carbon nanotubes [20,21,22], 
2d metallic [5,6,23,24,25,26,27] and semiconductor 
[28,29,30] films, various 3d disordered metals [26,27,31] 
and (Od) quantum dots [32,33,34,35,36]. Though di- 
mensions and parameters of these systems are differ- 
ent, the low temperature saturation of r v remains the 
common feature of all these observations. 

Is this ubiquitous saturation of t v an intrinsic or ex- 
trinsic effect? If intrinsic, decoherence of electrons at 
T = would be a fundamentally important conclusion 
which would shed a new light on the physical nature 
of the ground state of disordered conductors as well as 
on their low temperature transport properties. While 
extrinsic saturation of t v could be caused by a vari- 
ety of reasons, the choice of intrinsic dephasing mech- 
anisms is, in fact, much more restricted. There exists, 
however, at least one mechanism, electron-electron in- 
teractions, which remains important down to lowest 
temperatures and may destroy quantum interference 
of electrons even at T = [4,37]. 

In a series of papers [37,38,39,40] we offered a theo- 
retical approach that allows to describe electron inter- 
ference effects in the presence of disorder and electron- 
electron interactions at any temperature including 
the most interesting limit T — > 0. This formalism ex- 
tends Chakravarty-Schmid description [3] of WL and 
generalizes Feynman-Vernon-Caldeira-Leggett path 
integral influence functional technique [41,42,43,44] 
to fermionic systems with disorder and interactions. 
With the aid of our approach we have evaluated WL 
correction to conductance and electron decoherence 
time in the limit T — > and demonstrated that low 
temperature saturation of t v can indeed be caused 
by electron-electron interactions. Our results allowed 
for a direct comparison with experiments and a good 
agreement between our theory and numerous experi- 



mental data for t v in the low temperature limit was 
found [37,38,45,46]. In particular, for quasi-ld wires 
with thicknesses exceeding the elastic electron mean 
free path I at T — > our theory predicts t v oc D 3 , 
where D is the diffusion coefficient. This scaling is 
indeed observed in experiments for not very strongly 
disordered wires typically with D > 10 cm 2 /s (see 
Sec. 6 for more details). 

On the other hand, for strongly disordered structures 
with smaller values of D this scaling is not anymore 
fulfilled and, moreover, an opposite trend is observed: 
t v was found to increase with decreasing D [27,31,47]. 
This trend is not described by our expressions for t v 
[37,38]. Another interesting scaling was observed in 
quantum dots: saturated values t v were argued [36] to 
scale with the dot dwell times td as t v ~ td ■ Our the- 
ory [37,38] cannot be directly used in order to explain 
the latter scaling either. 

In order to attempt to reconcile all these observa- 
tions within one approach it is necessary to develop a 
unified theoretical description which would cover es- 
sentially all types of disordered conductors. It is worth 
pointing out that the technique [38,39,40] is formally 
an exact procedure which should cover all situations. 
However, for some structures, such as, e.g., quantum 
dots and granular metals, it can be rather difficult to 
directly evaluate the WL correction within this tech- 
nique for the following reasons. 

First of all, our description in terms of quasiclassi- 
cal electron trajectories may become insufficient in the 
above cases, and electron scattering on disorder should 
be treated on more general footing. Another - purely 
technical - point is averaging over disorder. In our ap- 
proach [37,38,39,40] it is convenient to postpone dis- 
order averaging until the last stage of the calculation. 
In some cases - like ones studied below - it might be, 
on the contrary, more appropriate to perform disorder 
averaging already in the beginning of the whole con- 
sideration. In addition, it is desirable to deal with the 
model which would embrace various types of conduc- 
tors with well defined properties both in the long and 
short wavelength limits. This feature will help to con- 
struct a fully self-contained theory free of any diver- 
gencies and additional cutoff parameters. 

Recently [48] we made a first step towards this uni- 
fied theory. Namely, we adopted a model for a disor- 
dered conductor consisting of an array of (metallic) 
quantum dots connected via junctions (scatterers) with 





Fig. 1 . Single quantum dot and a pair of time-reversed elec- 
tron paths. Fluctuating voltages Vi, and Vr are assumed 
to drop only across the barriers and not inside the dot. 

arbitrary transmission distribution of their conducting 
channels. This model allows to easily crossover between 
the limits of a granular metal and that with point-like 
impurities and to treat spatially restricted and spa- 
tially extended conductors within the same theoretical 
framework, as desired. Within this model in Ref. [48] 
we analyzed WL corrections to conductance merely for 
non-interacting electrons and included interaction ef- 
fects by introducing the electron dephasing time r v just 
as a phenomenological parameter. Systematic analysis 
of the effect of electron-electron interactions on weak 
localization within this formalism will be developed in 
this paper. This approach will allow to microscopically 
evaluate t v for all the structures under consideration. 
The structure of our paper is as follows. In Sec. 2 we 
will discuss qualitative arguments illustrating the role 
of scattering and interactions in electron dephasing. In 
Sec. 3 we introduce our model of an array of quantum 
dots and outline a general theoretical framework which 
is then employed in Sec. 4 and 5 for rigorous calcula- 
tions of the WL correction to conductance and electron 
decoherence time in the presence of electron-electron 
interactions. A detailed comparison of our results with 
numerous experiments performed in various disordered 
conductors is carried out in Sec. 6. A brief summary of 
our main results and conclusions is presented in Sec. 7. 



2. Qualitative arguments 

Before turning to a detailed calculation it is instruc- 
tive to discuss a simple qualitative picture demonstrat- 
ing under which conditions electron dephasing by in- 
teraction is expected to occur. 

Consider first the simplest system of two scatterers 
separated by a cavity (quantum dot, Fig. 1) The WL 
correction to conductance of a disordered system Gwl 
is known to arise from interference of pairs of time- 
reversed electron paths [3]. In the absence of interac- 



Fig. 2. Two quantum dots and a typical electron path. 
Fluctuating voltages Vl, V and Vr are again assumed to 
drop only across the barriers. 

tions for a single quantum dot of Fig. 1 this correction is 
evaluated in a general form [49] . The effect of electron- 
electron interactions can be described in terms of fluc- 
tuating voltages. Let us assume that the voltage can 
drop only across the barriers and consider two time- 
reversed electron paths which cross the left barrier 
(with fluctuating voltage VL(t)) twice at times ti and 
tf as shown in Fig. 1. It is easy to see that the voltage- 
dependent random phase factor exp(i J '/ VL{t)dt) ac- 
quired by the electron wave function \P along any path 
turns out to be exactly the same as that for its time- 
reversed counterpart. Hence, in the product $$* these 
random phases cancel each other and quantum co- 
herence of electrons remains fully preserved. This im- 
plies that for the system of Fig. la fluctuating voltages 
(which can mediate electron-electron interactions) do 
not cause any dephasing. 

This qualitative conclusion can be verified by means 
of more rigorous considerations. For instance, it was 
demonstrated [50] that the scattering matrix of the 
system remains unitary in the presence of electron- 
electron interactions, which implies that the only ef- 
fect of such interactions is transmission renormaliza- 
tion but not electron decoherence. In Ref. [51] a similar 
conclusion was reached by directly evaluating the WL 
correction to the system conductance. Thus, for the 
system of two scatterers of Fig. la electron-electron in- 
teractions can only yield energy dependent (logarith- 
mic at sufficiently low energies) renormalization of the 
dot channel transmissions [52,53] but not electron de- 
phasing. 

Let us now add one more scatterer and consider 
the system of two quantum dots depicted in Fig. 2. 
We again assume that fluctuating voltages are concen- 
trated at the barriers and not inside the cavities. The 
phase factor accumulated along the path (see Fig. 2) 
which crosses the central barrier twice (at times ti and 
t > ti) and returns to the initial point (at a time tf) 



o i[ v (t t )- v (t)} 



where <p/e = V(t) is the fluctuating 



voltage across the central barrier. Similarly, the phase 
factor picked up along the time-reversed path reads 
e iMt f +U-t)-v(.tf)], Hence, the overall phase factor ac- 
quired by the product *I>*I>* for a pair of time-reversed 
paths is exp(i<E> t ot), where 

$tot(ti, t f ,t) = tp(ti) - (f(t) - <p + (tf + U-t) + <p(t f ). 

Averaging over phase fluctuations, which for simplicity 
are assumed Gaussian, we obtain 



»*to* (*<,*/, t) 



■ I (*Lt(«i. */.*)} 



2F(t-t i )~2F(tf~t) + F(tf-t i )+F(tf+t i -2t) 



(1) 



where we defined the phase correlation function 

F(t) = ((#^(0)) 2 }/2. (2) 

Should this function grow with time the electron 
phase coherence decays and, hence, Gwl has to be 
suppressed below its non-interacting value due to 
interaction-induced electron decoherence. 

The above arguments are, of course, not specific to 
systems with three barriers only. They can also be 
applied to any system with larger number of scatter- 
ers, i.e. virtually to any disordered conductor where - 
exactly for the same reasons - one also expects non- 
vanishing interaction-induced electron decoherence at 
any temperature including T — 0. In the next sec- 
tions we will develop a quantitative theory which will 
confirm and extend our qualitative physical picture. 
We are going to give a complete quantum mechanical 
analysis of the problem which fully accounts for Fermi 
statistics of electrons and treats electron-electron inter- 
actions in terms quantum fields produced internally by 
fluctuating electrons. Below we will non-perturbatively 
evaluate WL correction Gwl for arrays of metallic 
quantum dots in the presence of electron-electron in- 
teractions which will be shown to reduce phase coher- 
ence of electrons at any temperature down to T = 0. 



3. The model and basic formalism 

Let us consider a Id array of N — 1 quantum dots 
connected in series chaotic quantum dots (Fig. 3) . Each 
quantum dot is characterized by its own mean level 
spacing S n . Adjacent quantum dots are connected to 
each other via barriers which can scatter electrons. The 
first and the last dot are also connected to the leads 




Fig. 3. Schematics of Id array of TV — 1 quantum dots 
connected via TV barriers. 

(with total resistance Rs), i.e. altogether we have TV 
scatterers in our system. Each such scatterer is de- 
scribed by a set of transmissions of its conducting chan- 
nels T^ (here k labels the channels and n labels the 
scatterers). Below we will focus our attention on the 
case of metallic quantum dots with the level spacing 
S n being the lowest energy parameter in the problem. 
The system of Fig. 3 will be described by the Hamil- 
tonian 



n,m= 1 

N-l N 



+ E^ d + E A 

n— 1 n— 1 



T 
n i 



(3) 



where C„ m is the capacitance matrix of the array, V n 
is the operator of electric potential on n— th quantum 
dot, 






tr \-L,R 

eV L ,n)a' h 



are the Hamiltonians of the left and right leads, Vl,r 
are the electric potentials of the leads fixed by the ex- 
ternal voltage source, 

Hf = J^ & i, a ,k( H n%i ~ eV n )a n ,a,i 

kl,a=-\,i 

is the Hamiltonian of n-th quantum dot and H n M is a 
randonr matrix which belongs to the orthogonal ensem- 
ble. Finally, electron transfer between adjacent n— 1-th 
and n-th quantum dots will be described by the Hamil- 
tonian 

hI= E /<*"*• M«-)*U-i(0*«,»(0+cc.]. 

Here the integration runs over the junction area. 

Note that in Ref. [48] we have already applied the 
model of Fig. 3 in order to analyze WL effects in the 
absence of electron-electron interactions. In that pa- 
per we have used the scattering matrix formalism com- 



bined with the non-linear <r-model. In order to incorpo- 
rate interaction effects into our consideration it will be 
convenient for us to describe inter-dot electron transfer 
within the tunneling Hamiltonian approach, as speci- 
fied above. We would like to emphasize that this choice 
is only a matter of technical convenience, since the ap- 
proaches based on the tunneling Hamiltonian and on 
the scattering matrix are fully equivalent to each other. 
For the sake of completeness, let us briefly remind 
the reader the relation between these two approaches. 
Consider, e.g., the n-th barrier between the n — 1-th 
and n-th quantum dots and define the matrix elements 
tim — (l\H„\m) between the I— th wave function in 
the n — 1-th dot and m— th wave function in the n- 
th dot. Electron transfer across this barrier can then 
be described by the set of eigenvalues of this matrix 
tu- These eigenvalues are related to the barrier channel 
transmissions T^ as ( see, e.g., [49]) 



T, 



(n) _ 47T \t k \ /8 n -iS„ 

' (l + 7r 2 |4| 2 /(5„_i<5„ 



(4) 



This equation allows to keep track of the relation be- 
tween two approaches at every stage of the calculation. 
To proceed we will make use of the path integral 
Keldysh technique. A brief sketch of this approach is 
outlined below. The time evolution of the density ma- 
trix of our system is described by the standard equa- 
tion 



~ / ,\ — iHt ~ iH t 

p(t) —e poe , 



(3) 



where H is given by Eq. (3) . Let us express the opera- 
tors e~ lHt and e lHt via path integrals over the fluctuat- 
ing electric potentials V F ' B defined respectively on the 
forward and backward parts of the Keldysh contour: 



. " = / DV n F T exp \ -i I dt'H [Vf{t'j\ 



t'H[\ 



e"" = J DV» T exp I i J dt'H \V» (i')J J • (6) 

Here T exp (T exp) stands for the time ordered (anti- 
ordered) exponent and the Hamiltonians H \V F (£')] , 
H \V B (£')] are obtained from the original Hamiltonian 
(3) if one replaces the operators V„(t) respectively by 
the fluctuating voltages V F (t') and V B (t'). 

The central object of our analysis is the effective 
action defined as 



iS[V F , V B ]=ln tr 



T exp 



x p T exp li I dt 



'h [v n F (f ; 

'H [v n B (t') 



-i I dt 
It) 



(•7) 



Since the operators H[V F (t')], H [V B (t')] are 
quadratic in the electron creation and annihilation op- 
erators, it is possible to rewrite the action in the form 



iS = 2 tr In G" 



(8) 



where G _1 is the inverse Keldysh Green function of 
the system. The operator G" 1 has the following matrix 
structure 



(Gl. 



G- x = 



T-t 



Ti 
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f 2 
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Gt 


f 3 







1 N-l 


f'- 1 


f N 








1 N 


Gr 



\ 



(9) 



V 

Here each quantum dot as well as each of the two leads 
is represented by a 2 x 2 diagonal block 



G n — 



idt - Hf + eV r 



I.) 



-idt + HT - eV, 



(10) 



while barriers are described by off-diagonal blocks T n , 
which have the form 



T n = 




(11) 



Below we will employ these general expressions in spe- 
cific situations of two quantum dots and of an array 
with a large number of quantum dots. The latter sys- 
tem will also serve as a model for spatially extended 
diffusive conductors. 



4. System of three barriers 

Let us first consider the system of three scatterers 
iV = 3 (two quantum dots) . This system is important 
both in the context of electron dephasing in quantum 



dots [32,33,34,35,36] and as a relatively simple exam- 
ple which illustrates all significant features of the effect 
of electron-electron interactions on WL in more com- 
plicated systems. 

Here we will treat interaction effects practically with- 
out any approximations. For the sake of simplicity we 
will assume that the central (i.e. second) barrier is a 
tunnel junction with Tk <JC 1 for all its transmission 
channels. This assumption allows one to expand the 
action (8) in powers of the parameter t 2 (r). We will 
also assume that dimensionless conductance of the cen- 
tral barrier g 2 = 2-k /e 2 R 2 is much lower than those of 
the first and third barriers which, in turn, strongly ex- 
ceed unity, i.e, gi,gs 3> g 2 , 1. The effective expansion 
parameter in this case is </f/<?i33 *C 1- Then the WL 
correction to the conductance in the absence of inter- 
actions takes a particularly simple form [48]. E.g., for 
fully open outer barriers we have 



.^ array 



-.(()) 



2e 2 



1 + M + r (0) 



<)i 



gl 



i" gi9i 



(12) 
(13) 



The first term in Eq. (12) defines the contribution of the 
first and third scatterers, while the term G WL comes 
from the central junction. Experimentally one can ac- 
cess G]J L attaching additional voltage probes to the 
quantum dots. The aim of our subsequent analysis is 
to demonstrate how the result (13) is modified in the 
presence of electron-electron interactions. 



4.1. WL correction in the presence of interactions 

Setting N = 3 in Eq. (3) and defining the fluctuating 
phases across the barriers <p^' B = /. dt' [eV^' B (t') — 
eV„li (t')] we perform the gauge transformation which 
yields the new expressions for the blocks G„ and T n : 



' ' ■ 



T n 



idt - H% 




—t n (r)e 




qd 



-id t + Hf 



t„(r)e 



(14) 



(15) 



Since the central barrier transmission is small we can 
expand the action (8) in powers of the parameter £ 2 (r). 



Proceeding with this expansion up to the fourth order 
we get 

iS = iSs[<pf£] + iS { 2 2) + iS { 2 4) . (16) 

Here 

iSs&if] = iSc + *S cxt + 2tr mlG- 1 ])^^ (17) 

is the effective action at zero transmission of the 
second barrier which also includes terms describing 
capacitances (iSc) and the external circuit (iS ex t), 
the term iS 2 oc t 2 (r) is the standard Ambegaokar- 
Eckern-Schon (AES) action [43] and iS ( 2 ] oc ti(r) 
contains information which will allow us to evaluate 
the WL correction to the system conductance. The 
corresponding expression reads 



iS i 4) =- Yl [ dti...dU f 



i,3,k,l=F,B 



dtA I dxi . 



. dxA 



xG 1 , lJ (A 1 ;X 2 )(-l) 3 



-ifi(t 2 ) 



t(x 2 ) 



xG 2 , ]k (X 2 ;X :i )(-l) k e^ (t3) t(x 3 ) 
x G 1M {X 3 - X 4 )(-l) l e-'^ M t(xi) 
xft ! i,(l4;Ii)(-l)V^ |t ''t( I1 ). 



(18) 



Here we use the convention ( — 1) F = — 1, (— 1) B — 1, 
X = (t,x), and G r are 2x2 matrix Green-Keldysh 
functions in the first and the second quantum dots (r — 
1,2): 

/ (T^(x 1 )^+(x 2 )> -(4,i(x 2 )MXi)) . m . 

lG r = I | (19) 

<</v(Xi)$(X 2 )> {fi, r (X 1 )^ r (X 2 )} 

Here we will set the Green functions Gi and G 2 equal 
to their equilibrium values G r = G^F\ — F 2 G r , where 
G B ' are retarded and advanced Green functions, 



F, = 



h(E) -f(E) 
ME) -f(E) i 



F 2 = 



f(E) f(E) 
-h(E) -h{E) 



(20) 

f{E) is the Fermi function and h(E) = 1 - f(E). This 
choice is sufficiently accurate for the problem in ques- 
tion. We will return to this point towards the end of 
this section. 

Our next step amounts to averaging the products of 
retarded and advanced propagators in the action (18) 
over disorder in each of the two dots separately. This 
averaging can be accomplished, e.g., with the aid of 



the non-linear cr-model or by other means. For the first 
(left) dot we obtain (cf., e.g., [54]) 



{G H (X 1 ,X 2 )G?(X 3 ,X 4 )} = 
27rJViViw(|ri — r4|)w(|r 2 — r 3 \) 

x Di f ii - t 2 ; 5 ' o J ^ ~ ~ ^ 

+ 27riViViw(|ri — r3|)w(|r2 — r4|) 

x Ci Ui -tr, — - — , — - — <J(*i -t 2 + t 3 -U) 



+(G?(X 1 ,X 2 )}{G?(X a ,X 4 )}, 



(21) 



where Vi is the volume of the first (left) quantum dot, 
Ni, D\(t,r,r') and Ci(t,r,r') are respectively the 
density of states, the diffuson and the Cooperon in the 



first (left) dot, w(r 



-rill 



sin kpr/kFr, Uf and I are 



respectively the Fermi wave vector and elastic mean 
free path. The same averaging procedure applies to the 
second (right) dot. In the absence of magnetic field 
both the diffuson and the Cooperon satisfy the same 
diffusion equation 



m 



DV 2 ) Di (t, r,r') = Vi6(t)6{r - r 



(22) 



with the appropriate boundary conditions at the con- 
tacts. 

Averaging the action (18) over disorder we will col- 
lect only the terms proportional to the product of the 
Cooperons C\ and C 2 and ignore all other contribu- 
tions which are unimportant for the problem in ques- 
tion. In this way we arrive at the action Swl which 
describes weak localization effects in our system. Ig- 
noring for simplicity the coordinate dependence of the 
Cooperons G\ , C 2 inside the dots, we obtain 



iSwl 



. g 2 SiS 2 

4-k 2 



dti . . . dti / dr\dr 2 Ci(ti — n) 



xC 2 (fe-r 2 )e i ^ + ^)-^ + fe)+^ + ('3)-^(*4)] sin ^i*ll 



h(n - t 2 )e 



= + /( Tl _ t2 y 2 



h{r 2 - t 3 )e i ^ 3 ~f{ti + t 3 -U- n) 
/ (r 2 - t 3 )e~ 1 ^^ h(h + t 3 -U- Tl) 



i V (*4) „. , -je_J*4) 

e 2 /(-ti + £ 2 + £4 - r 2 ) + e 2 



x ft(-ti +i 2 +*4 -r 2 ) 



+ {1^2, ¥ j ± -h--^}. (23) 



Here we defined mean level spacing 5i, 2 for both 
dots, introduced "classical" p + = (ip F + <p B )/2 and 
"quantum" p~ — ip F — ip B phases and made use of 
the Fourier transforms of the Fermi function f(t) = 
J(dE/2Tv)f(E)e- iEt = 8{t) - h(t). The resulting 
action (23) fully accounts for the effects of electron- 
electron interactions on WL via the fluctuating phases 
„±. 

In order to find the WL correction to the current 
across the central barrier Iwl we employ the following 
general formula 



Iwl 



_2 ± SiSwL 

D p — 

dtp 



[^ 



iSs+iSj 



(24) 



The task at hand is to combine Eqs. (23) and (24) 
and to average over the phases tp by evaluating the 
path integral in (24). The contributions Sc and S^xt 
in (17) are quadratic in the fluctuating phases pro- 
vided an external circuit consists of linear elements. 
The remaining contribution to Ss (17) which describes 

transfer of electrons across the first and third barri- 

(2) 
ers as well as AES action of the second barrier S 2 

are in general non-Gaussian. However, in the interest- 
ing limit gi t s 3> I phase fluctuations can be consid- 
ered small down to exponentially low energies [55,56] 
in which case it suffices to expand both contributions 
tr ln[G _1 ]| . and S 2 up to the second order ip . 

Furthermore, Gaussian approximation for the first of 
these terms becomes essentially exact in the limit of 
fully open outer barriers with gi t3 3> 1 [57]. 

We conclude that the integral (24) remains Gaussian 
at all relevant energies and can easily be performed. Af- 
ter straightforward algebra we arrive at the final result 



Re / dEdu)idu) 2 du> 3 C 2 (— u> 2 ) 



T eg 2 Si5 2 

x Ci(-u> 3 )h{E - u> 2 )f(E + eV + uj 3 - wi) 
x [f{E + eV-u>i)h{E)P 1 (wi,W2,w 3 ) 

+ /(£ + eV - Wi)f(E)P 2 (u>i, W2,W 3 ) 

+ /i(£ + eV - wi)ft(S)P 2 (wi , ^3,^2) 
+ fc(E + eV - u>i)f{E)P 3 (u>i, u> 2 , w 3 )] 
-{^->-V}. (25) 

Here V is the average voltage across the central bar- 
rier, Ci,2 are the Fourier transforms of the Cooperons 
Ci,2(t) and the functions Pj (j = 1, 2, 3) are defined as 

P i(wi ) = I df ff t3 e^+^+^'slp.^ (26) 

where 
Pj(ti,t 2 ,t 3 ) =exp[-r{ti,t2,t3)}Qj{ti,t2,t 3 ), (27) 

^(ti, t 2l is) = P(ti + ts) + F(t 3 ) + F(ti + t 2 ) 
+F(t 2 ) - F{ti +t 2 + t 3 ) - F(t 2 - t 3 ), (28) 

and Fit) = {(ip(t) — <p(0)) 2 }/2 coincides with the phase 
correlation function (2). The terms Qj read 

Q 1 _ e -»[^(*2) + ^(«3) + ^(l*2-t 3 l)] 

x {2e ^[K(|tl+t2+t3l)+i:f(tl+t3)+K(tl+t2), 

_ e »[A-(ti+t 2 +*3)+«'(l*l+*3l)+-K'(l*l+*2l)] I 
i[K(\t 1 +t 2 +t 3 \)-K(t 2 )-K(\t 3 \)} 



, JlK(t 1 +t 3 )-K(\t 1 +t2\)-K(t3-t 2 )] 

"■ e ) 

) g = e l[-K(tl+t2+*3)--K(l*2l)-Jf(|t3l)] 
, -<[^(tl+«3)+^(*l+t2)-if(l*3-*2l)] 



(29) 



where -K"(t) — i( [0(0), >p{t)] ) is the response function. 
Eqs. (25)- (29) fully determine WL correction to the 
current in our system. The non-interacting result [48] is 
reproduced by the first line of Eq. (25) , while its second 
and third lines exactly account for the effect of inter- 
actions. We observe that the whole effect of electron- 
electron interactions is encoded in two different corre- 
lators of fluctuating phases F(t) and K(t) defined as 

(t) = / 27 W C ° 2T Re [ (W) ] ^ . ( 3 °) 



where Z(ui) is an effective impedance "seen" by the 
central barrier. 

Let us recall that these correlation functions play a 
central in the so-called P(iJ)-theory of Coulomb block- 
ade in tunnel barriers [43,58]. It turns out that exactly 
the same correlators also describe the effect of electron- 
electron interactions on weak localization. Below wc 
will demonstrate that the phase correlation function 
(30) is responsible for electron dephasing while the re- 
sponse function (31) describes the Coulomb blockade 
correction to WL. 

4.2. Dephasing 

Let us now turn to the analysis of the above general 
results. To begin with, we notice that both functions 
(30) and (31) are purely real and, hence, \Qj\ < 1. Fur- 
thermore, at sufficiently long times |t| > trc (where 
trc is an effective PC-time of our system to be defined 
later) we obtain 



2 / sintiTrTi 
gz V kTtrc 
K(t)~ — sign*, 



(32) 
(33) 



where gz = 2ir/e 2 Z(Q) = g + g 2 , g : = g 1 : + g 3 i + 
e Rs/2ty and 7 ~ 0.577 is the Euler constant. The 
correlation function F(t) grows with time [59] at any 
temperature including T = 0. In contrast, the response 
function K (t) always remains small in the limit gz ^> 1 
considered here. Hence, the combination (28) should be 
fully kept in the exponent of (27) while the correlator 
K(t) can be safely ignored in the leading order in 1/ gz- 
Then all Qj = I, the Fermi function f(E) drops out 
from the result and we get Iwl — GwlV , where 



Ch 



where 



e 2 g 2 t S 1 S 2 

8tt 3 
x C\(t 2 )C 2 (t 3 ) 



dt 2 dt 3 e- T(t2M) 



(34) 



. ;' du) _ r™ si siliiut 



(31) 



F{t 2 , t 3 ) = 2F(t 2 ) + 2F(t 3 ) - F{t 2 - t 3 ) - F{t 2 + t 3 ). 

(35) 
Identifying t 2 — tf—t and t 3 — t — i, in Eq. (35) we ob- 
serve that the exponent exp(— J-(t 2 , t 3 )) exactly coin- 
cides with the expression (I) derived from simple con- 
siderations involving electrons propagating along time- 
reversed paths in an external fluctuating field. Thus, 



we arrive at an important conclusion: In the leading 
order in 1/gz the WL correction Gwl is affected by 
electron-electron interactions via dephasing produced 
only by the "classical" component ip + of the fluctuat- 
ing field which mediates such interactions. This effect 
is described only by the phase correlation function F(t) 
(30). At the same time, fluctuations of the "quantum" 
field ip~ turn out to be irrelevant for dephasing and 
may only cause a (weak) Coulomb blockade correction 
described by the response function K(t) (31). This lat- 
ter effect will be analyzed in the next subsection. 

In order to simplify our consideration let us assume 
that two quantum dots are identical, i.e. gi = g-a = 
j> 1,32 and Si — 62 — S, and set Rs — > 0. In this 
case the Cooperons are C\(t) = C 2 (t) = e -*/ TD -*/ r » , 
where td — An/gS is the dwell time for each of the two 
dots and th oc 1/H 2 is the dephasing time due to the 
external magnetic field H which can be applied to the 
system. An effective impedance Z(cu) takes the form 



Re Z(lu) = — 
■Kg 



tyS(uj) 



1+UJ 2 T 2 



TD + TRC 



(36) 



where 1/r = 1/td + 1/trc, trc = Tv/gEc, Ec — 
e 2 /2(C* + C 3 + 2Cj) and C, C.j and C 3 are the ca- 
pacitances of respectively left(right) barriers, the cen- 
tral junction and the gate electrode. Substituting the 
Cooperons Ci^if) into Eq. (37) we arrive at the final 
expression for the WL correction Gwl (T) in the pres- 
ence of electron-electron interactions 



Gwl 



e g 



8tt 3 



(t 2 +t 3 )( Tp +T H ) 

dt2dtze t d^h 



(37) 



With a good accuracy the double time integral in Eq. 
(37) can be replaced by a single one, in which case the 
magnetoconductance Gwl can be expressed in a much 
simpler form 



G\ 



Til 



G 



(o) 



Td(th + To) 



dte- 2F(t/2) e 



t/rn-t/T 



'(SB) 



The result (37) for H = is plotted in Fig. 4. We 
observe that electron-electron interactions always sup- 
press Gwl{T) below its non-interacting value (13). 
This is a direct consequence of interaction-induced 
electron dephasing described by the correlation func- 
tion F(t). 
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Fig. 4. Temperature dependence of WL correction Gwl 
in the presence of electron-electron interactions for 
td/trc = 100. 

Let us define u — td/trc = 4Ec/S and consider the 
limit of metallic dots u ^> 1. At g/TD < T < 1/trc 
and for j>8we find 

(2tt/m) 8/9 



IG 



(") 



- 1 1 - 2 ^ 



(•kTtd) 1 - 8 ^' 



(39) 



whereas at Ttd 5; 1 the WL correction saturates to 



Gwl/\G wl \ 
Gwl/\G wl \ 



-(2/n) 8/9 , g>8, 
-S/2M, 1<5<8, 



(40) 



Let us compare the magnitude of the WL cor- 
rection (40) to that of the leading contribution to 
the system conductance evaluated in the presence of 
electron-electron interactions. This contribution reads 

[60,43,58] 

d/(T = 0) 1 [e-"eVTRcY h 

dv r 2 r (1 + 4/0) ' v ' 

where Y(x) is the gamma-function. Setting eV ~ 1/td 
in the above expression, we observe that the WL cor- 
rection (40) remains much smaller than Eq. (41) by the 
parameter ~ g2/gig3U 4 ' 9 . 

Let us phenomenologically define the electron deco- 
herence time t v by taking the Cooperons (for H — + 0) 
in the form Ci, 2 (t) = e" t/TD " t/T ". This definition ob- 
viously yields [48] G W l/G$ l = (1 + t d /t v )- 2 . Re- 
solving this equation with respect to t v , we obtain 

-1 



?/td 



G W l/Gwl 



(42) 



Substituting the result (39) into Eq. (42) at sufficiently 
high temperatures g/TD < T < 1/trc and in the limit 
g» 1 we obtain 
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Fig. 5. Temperature dependence of the electron dephasing 
time r v in the presence of interactions for t^/trc = 100. 



t v {T) ~ ^/gr D 



Ec) 



4/9 



(Tro) 



4/S-1/2 



(43) 



Combining Eqs. (40) and (42) at lower temperatures 
Ttd < 1 and for j>8we find 

TD 



T<p(T) — T V () — 



(2E c /S) 4 /9 - 1 ' 



In the limit of large g 3> 1 this expression yields 
gr D tv 



(44) 



(45) 



4\n(2E c /S) 8\n(2E c /S)' 

The results for the electron decoherence time t v (T) 
(42) are also plotted in Fig. 5 for different values of g. 
We observe that at higher temperatures t v {T) shows a 
(/-dependent power law dependence on T and eventu- 
ally saturates to the value (45) at lower T. This is ex- 
actly the behavior observed in a number of experiments 
with quantum dots [32,33,34,35,36]. We will postpone 
a detailed comparison between our theory and experi- 
ments to Sec. 6. 



4.3. Perturbation theory and Coulomb correction to 
weak localization 

Let us now analyze the role of the response function 
K(t) which was disregarded in the previous subsection. 
For this purpose let us expand the general expression 
for the current (25) to the first order in the interaction, 
i.e. to the first order in both F(t) and K(t). We get 



Iwl — (j wl 



V + 6I^ L (V)+SI« L (V), 



(46) 



F e g 2 5iS 2 V 



duo HeZ(uj) , lu 

— L coth — - 

2tt u) 2T 



8-K 3 J 2lT LU 

x [2C , i(0)(5 2 (w) + 2C , i(w)C , 2 (0) - 2C 1 (0)C 2 (0) 
-CiHCaH-Cit-o/lC'aH], (47) 



„k e g 2 5 1 8 2 

OlwL = 



dm W(u;,V) 
16tt 3 J 2tt lo 

x [2d(Q)C 2 (uj) + 2C 1 {lu)C 2 {0) - 
+ iImZ(uj)Ci(uj)C 2 (uj)}. 



{ReZ(oj) 



(48) 



(o> — eV) coth ■ 



eV 



where 



Here we defined the function 

W = (u> + eV) coth - + eV - v 

The above expressions allow to make several impor- 
tant observations. To begin with, we note that the term 
SIwl is linear in the bias voltage while &Iwl is non- 
linear in V . The physical meaning of these two terms is 
entirely different. While, as we already know, the cor- 
rection &I F r L describes electron dephasing, the term 
SlyvL i s nothing but the Coulomb blockade correction 
to Gwl (37). The latter conclusion is supported, for in- 
stance, by the observation that at large voltages SI^l 
tends to a constant offset value, which is a typical sign 
of Coulomb blockade. 

Just for the sake of illustration, let us for a moment 
assume that the environment remains Gaussian at all 
energies. In this case dI^r L /dV diverges at V, T — > 
indicating the importance of higher order perturba- 
tive terms in the low energy limit. Making use of an- 
alytical properties of the functions Pj(u)\,u) 2 ,uj?,) one 
can exactly sum up diverging perturbative series to 
all orders and quite generally demonstrate that in the 
limit of zero voltage and temperature the WL correc- 
tion tends to zero, dlwL/dV — ► 0, implying complete 
Coulomb suppression of weak localization in this limit. 
This strong Coulomb blockade of WL can be recov- 
ered only if one fully accounts for the response func- 
tion K(t). We also note that at large conductances 
weak Coulomb blockade may turn into strong one only 
at exponentially small energies [55,56] typically well 
below the inverse electron dwell time 1/td, otherwise 
Coulomb suppression of WL remains weak. 

Turning back to the first order terms (47) and (48), 
we observe that in the linear response regime V — > al- 
most all contributions from 6 Iwl an< i &Iwl cancel each 
other exactly in the limit T — > 0. This (partial) can- 
cellation of the so-called "coth" (47) and "tanh" (48) 
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terms is a general feature of the first order perturbation 
theory. For instance, in diffusive conductors it has been 
observed by various authors [61 ,62]. This cancellation 
is sometimes interpreted as a "proof" of zero dephas- 
ing of interacting electrons at T = assuming that the 
same "coth-tanh" combination should occur in every 
order of the perturbation theory [63]. Our exact result 
(25) demonstrates, that this assumption is incorrect. 
Partial cancellation between dephasing and Coulomb 
blockade terms occurs only in the first order and is of 
little importance for the issue of electron dephasing. 
Actually, the combination 1 - 2f(E) = tanh(E/2T) 
enters only in the first order, while the exact expres- 
sion (25) depends on f(E) and 1 — f(E) and does not 
contain the combination 1 — 2f(E) at all. This obser- 
vation is fully consistent with our general result for the 
WL correction [40] expressed in terms of the matrix 
elements of the operators p and 1 — p, where p is the 
electron density matrix. In fact, the analysis presented 
above is a way to explicitly evaluate these matrix ele- 
ments for a particular case of two quantum dots. 

Finally, let us estimate the corrections SI W l an d 
SI WL . F° r niax(eV^T) > 1/td we obtain 



dI WL Se 2 g 2 2 
dV irg 3 


7r7Y D + 21n-^-l 




(49) 


dI WL 

dV ~ 


f 20 
3 
6 

. "9 3 


g 6 sZ> eV « T > 

;Zr eV>>T > 


(50) 



i.e. at Ttd ~ 1 the ratio of the terms (50) and (49) 
is of order ~ 1/ \n(AEc / 5) <C 1 and it decreases fur- 
ther with increasing temperature as oc 1/T 2 t d . This 
estimate demonstrates that dI WL /dV can be ignored 
as compared to the main contribution dI^ VL /dV . This 
conclusion remains valid also beyond the first order 
perturbation theory. Indeed, keeping the function F(t) 
in the exponent of Eq. (25) and expanding this result 
only in powers of K(t), at Ttd > 1 one finds 



4.4. Summary of approximations 

For clarity, let us again summarize all the approxi- 
mations used in the above analysis: 

(i) Throughout our calculation we have used the 
equilibrium form of the distribution function matrices 
-Fi,2 (20) which effectively implies neglecting the de- 
pendence of the Green functions 61,2 on the phases <pi 
and y>3. This is accurate except at energies well below 
the inverse dwell time 1/tzj. The WL correction (34)- 
(35) saturates at energies ~ 1/td and, hence, is totally 
insensitive to this approximation. The Coulomb term 
dlwL/dV should be treated somewhat more carefully 
in the limit T,eV < 1/td- For Rs — > this treat- 
ment also yields saturation of the Coulomb correction 
to WL at energies of order 1/td, exactly as in the 
case of the interaction correction to the Drude conduc- 
tance [50,64]. In other words, at all T, eV < 1/td the 
term dI^ L /dV remains constant of order ~ Gwh/g <C 
Gwl- Thus, the above approximation is completely 
unimportant for any of our conclusions. 

(ii) We have assumed 31,3 S> 1 and 31,3 3> 32- The 
first inequality just implies that our structure is metal- 
lic, while the second one is only a matter of technical 
convenience. Obviously, our results for t v remain qual- 
itatively the same also for #1,3 ~ g 2 . 

(Hi) For completeness, let us also mention that in the 
final expressions for Gwl, e.g., in Eqs. (39) and (40), 
we have disregarded small renormalization of the con- 
ductance g due to Coulomb effects [52,53]. This renor- 
malization is strictly zero for fully open outer barriers, 
otherwise it can be trivially included into our final re- 
sults. 

Thus, our treatment presented in this section does 
not contain any approximations which could influence 
our main results and conclusions. In the limit (ii) our 
general results for the WL correction to the I — V curve 
in the presence of electron-electron interactions (25)- 
(29) are exact. 



dT K 

al WL 

dV 



(7 y 



gTT D 



<Gv 



(51) 



5. Arrays of quantum dots and diffusive 
conductors 



thus supporting our conclusion about relative unim- 
portance of the Coulomb blockade correction to Gwl 
for the problem in question. 



One of the main conclusions reached in the previous 
section is that the electron decoherence time is fully 
determined by fluctuations of the phase fields ip + (and 
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the correlation function F(t)), whereas the phases <p~ 
(and the response function K(t)) are irrelevant for t v 
causing only a weak Coulomb correction to Gwl ■ This 
conclusion is general being independent of a number of 
scatterers in our system. Note that exactly the same 
conclusion was already reached in the case of diffu- 
sive metals by means of a different approach [38,39,40]. 
Thus, in order to evaluate the decoherence time for in- 
teracting electrons in arrays of quantum dots it is suf- 
ficient to account for the fluctuating fields V + totally 
ignoring the fields V~ . The corresponding calculation 
is presented below. We will specifically address the ex- 
ample of Id arrays and argue that the final result for the 
zero temperature decoherence time r^o actually does 
not depend on the dimensionality of the array. We will 
also demonstrate how our present results for t v match 
with those obtained previously for weakly disordered 
metals [37,38,391. 



This solution can be represented in the form C„m(t) = 
C^l(t) - C™t(t), where 



r* bulk (t\ 



N ^ J 2tt 

9=1 J 



■nq{n — m) 
N 



(55) 



-iu> + 



In the limit of large N the term Cjj"'^(t) can be 
safely ignored and we obtain C nm (t) ~ C^^it). Let 
us express the contribution C„t m (i) as a sum over the 
integer valued paths v(t), which start in the m— th dot 
and end in the n— th one (i.e. z/(0) = m, v{t) = n) 
jumping from one dot to another at times tj. This ex- 
pression can be recovered if one expands Eq. (55) in 
powers of r^, 1 cos[7rg/iV] with subsequent summation 
over q in every order of this expansion. Including addi- 
tional phase factors acquired by electrons in the pres- 
ence of the fluctuating fields Vj 1 " , we obtain 



5.1. Id structures 

Let us consider a ID array of N — 1 quantum dots 
by N identical barriers as shown in Fig. 3. For simplic- 
ity, we will stick to the case of identical barriers (with 
dimensionless conductance j> 1 and Fano factor f3) 
and identical quantum dots (with mean level spacing 
S and dwell time td — 2ii/5g). The WL correction to 
the system conductance has the form [48] 



Gwl = 



4tt 2 N 2 



E 



dt 



x{/3[C„_l,„(t)+Cn,n-l(*)] 

+ (l-f3)[C nn (t)+C n -i,n-i(t)]}. (52) 

The Cooperon C„m{t) is determined from a discrete 
version of the diffusion equation. For non-interacting 
electrons and in the absence of the magnetic field this 
equation reads 

OLyjim £(-^nm ^n— l,ra ^n+1,77 



+ 



dt 2 TD = S - 5{t) - (53) 

The boundary conditions for this equation are C nm = 
as long as the index n or m belongs to one of the bulk 
electrode. The solution of Eq. (53) with these boundary 
conditions can easily be obtained. We have 



2 N ~ 1 f 

C n,m(t) = -T 2_^ j 

" 0=1 J 



dui 
2m' 



sin sp sin Zjp 



-ioj + ■ 



(54) 



Cnm\t) — 2_^ / , 

fc— \n— m\ i/(t) 

i 






X 



(2r D )*J dtk J dtk ~ l 



d±2 \ dt\ 

il Jo 



t-'fc t k-*k-l t 2 -t 2 *2~*1 *1 

x e t d e t d . . . e t d e t d e t d 

x ^p[ijyr[eV+ T) (r)~eV+ t _ T) (r)]Y (56) 



Averaging over Gaussian fluctuations of voltages V + 
and utilizing the symmetry of the voltage correlator 
<V+(n)V+(7*)) = {V+WV+fa)), we get 

oc v{t)=n 

Cnm{t) = 2_^ Z-^l 

k=\n-m\ v{t) "(0)=m 



,-t/T D ft 



(2r D )* J dtk J dtk ~ l 



dti \ dt\ 

i) Jo 



x exp | - e 2 J* dn J* dr 2 [(V+ ti) (t 1 )V+ T2) (t 2 )} 
-<^n)W^)W)]| (57) 



The correlator of voltages can be obtained, e.g., with 
the aid of the cr-model approach employed in Ref. [48]. 
Integrating over Gaussian fluctuations of the Q-ficlds 
one arrives at the quadratic action for the fluctuating 
fields V + which has the form 
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2 ^ f dco gr D e 2 ( 1 



iV 



cos ^H ojcoth ^r 



2tt 



9=1' 

xV-(u)V-(-u). 
Here we defined 

JV- 
r± 



v?M = E 



U,V*+(l-C06#) 



disin^e^'V^i). 



(58) 



(59) 



The action (58) determines the expressions for both 
correlators (V+V+) (F-function) and (V+V-) (K- 
function) responsible respectively for decoherence and 
Coulomb blockade correction to WL. Since our aim is 
to describe electron decoherence, only the first out of 
these two correlation functions is of importance for us 
here. It reads 



<K + (ti)\/+(t 2 )) 



N 



JV-l 

E 



^- i 



9=1 

COS -rf ) sin 
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2n" 



-iw(t 1 -t 2 ) 



irqn 



TT </ III 



J4C(l-cos^)+C 9 + ^- 

T D U) coth ^r 

W 2 rB + (l-cos^) 2 ' 



— ilOTjj+l — COS -0- | 



(60) 



In the continuous limit N ~^> 1 and for sufficiently low 
frequencies w <JC 1/td both correlators (l /+ V r+ ) and 
(V + V~) defined by Eq. (58) reduce to those of a dif- 
fusive metal [38] . 

To proceed let us consider diffusive paths v(t), in 
which case one has 

1 N - 1 
(K^inW^ir^^j^-j E (^ + (n)\/+(r 2 )) 

XD nm {\Tl-T 2 \), (61) 

where D nm (r) is the diffuson. For H — *• it exactly 
coincides with the Cooperon for non-interacting elec- 
trons (54), D nm (t) = C*i°l(i), i.e. 



D nm (t) = - e / ^ 



9=1 



(L,' 



2tt 



7rqn ■ iram 

sin -^- sin -^ r 
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(62) 



Substituting Eq. (61) into (57), we obtain 

a m (f)«C<S,(t)e- rw , (63) 



where 

2 JV-l ,. t 

Ht) = jf—[ E / dUdt 2 {V+ {t 1 )V+{t 2 )) 

n,m=l , ' 

x [£>„ m (|ti - i 2 |) - £»„ m (|t - ti - *a|)] (64) 

is the function which controls the Cooperon decay in 
time, i.e. describes electron decoherence for our Id ar- 
ray of quantum dots. The WL correction Gwl in the 
presence of electron-electron interactions is recovered 
by substituting the result (63) into Eq. (52) . 

Note that in the limit of very large N and for t ^$> 
td Eq. (64) (combined with (60) and (62)) reduces to 
Eq. (22) of Ref. [39] derived for metallic conductors by 
means of a different technique. Replacing the sum in 
(64) by momentum integration, in the same limit one 
arrives at the result which matches with Eq. (23) of 
[39] provided one sets the charging energy terms equal 
to zero. 

Since the behavior of the latter formula was already 
analyzed in details earlier [39], there is no need to re- 
peat this analysis here. The dephasing time t v can be 
extracted from the equation J-(t v ) = 1. From Eq. (64) 
with a good accuracy we obtain 



N-l „ 

- E J dr{V+(r)V+(0))D nrn (r). (65) 



N- 



Combining this formula with Eqs. (60) and (62), in the 
most interesting limit T — > and for td ^> R(AC + C g ) 
we find 



1 1 N ~ 1 

Tipo ~ 2 9 t d {N -1) ^ n 



2c 2 



t v o 2gr D {N-l) 
which yields 



6 (AC (1 - cos #) + C g ) 



T V 



2(JT D 



An 



(66) 



ln{AE c /S) Sln(AEc/8) 

where E c = e 2 /2C a for C g > C and E c = e 2 /AC in 
the opposite case C g <§C C. 

We observe that apart from an unimportant numer- 
ical factor of order one the result for t v o (66) derived 
for Id array of quantum dots coincides with the ex- 
act result (45) derived in the previous section for the 
case of two quantum dots. Thus, we arrive at an im- 
portant conclusion: In the low temperature limit the 
electron decoherence time t v o is practically indepen- 
dent of the number of scatterers in the conductor (pro- 
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vided the latter number exceeds two) and is essentially 
determined by local properties of the system. 

In order to determine the dephasing length L v = 



'Dt v let us define the diffusion coefficient 



D = 



df_ 
2t d 



d 2 g5 

4tt ' 



(67) 



where d = V 1 ' 3 is the average dot size. Combining Eqs. 
(66) and (67), at T = we obtain 



L<po = s/Dt v0 = d\Jg/\n(4E c /5). 



(68) 



At non-zero T thermal fluctuations provide an addi- 
tional contribution to the dephasing rate 1/t v . Again 
substituting Eqs. (60) and (62) into (65), we get 

— + ^mm{N,N^}, (69) 

t v o 6g 



r v (T) 



where N v = L^/d ~ \Jt v /td is the number of quan- 
tum dots within the length L v . We observe that for suf- 
ficiently small N < N v (but still N ^$> 1) the dephasing 
rate increases linearly both with temperature and with 



the number N. At larger N > -J g/ ln[4Ec/8] and/or at 
high enough temperatures N v becomes smaller than N 
and Eq. (69) for r v should be resolved self-consistently. 
In this case we obtain 



(3gV^/^T) 2/3 , 



(70) 
(69) 



which matches with the AAK result [65]. Eq. 
also allows to estimate the temperature T* ~ 
3g/[TTTipo min{iV, N v }] at which the crossover to the 
temperature-independent regime (66) occurs. We find 









31n[4E c /S] 
2-kNtd 

3\n 3/2 [4E c /5] 

2-KTD^/g 



N< 



\n[AE c /5] ' 



N> 



\n[4E c /S] 



(71) 



5.2. Good metals and granular conductors 



The above analysis and conclusions can be gener- 
alized further to the case 2d and 3d structures. This 
generalization is absolutely straightforward (see, e.g. 
[48]) and therefore is not presented here. At T —> one 
again arrives at the same result for r^o (66). 

Now we discuss the relation between our present re- 
sults and those derived earlier for weakly disordered 
metals by means of a different approach [37,38,39]. Let 
us express the dot mean level spacing via the average 



dot size d as 8 — 1/Nod 3 (where No = mpF/2ir 2 is the 
electron density of states at the Fermi level). Then we 
obtain 

g 



Below we consider two different physical limits of (a) 
good metals and (b) strongly disordered (granular) con- 
ductors. For the model (a) we assume that quantum 
dots are in a good contact with each other. In this 
case g scales linearly with the contact area A = 7c? 2 , 
where 7 is a numerical factor of order (typically smaller 
than) one which particular value depends on geometry. 
For weakly disordered metals most conducting chan- 
nels in such contacts can be considered open. Hence, 
g — p%A/2it and 



D — 7«F<i/4, 



(73) 



i.e. D (x d. Comparing this estimate with the standard 
definition of D for a bulk diffusive conductor, D = 
vfI/3, we immediately observe that within our model 
the average dot size is comparable to the elastic mean 
free path, I ~ yd, as it should be for weakly disordered 
metals. 

Expressing r^o (66) via D, in this limit we get 



64 m 2 



D A 



where m is the electron mass and D c \ is constant which 
depends on Ec- Estimating, e.g., Ec ~ e 2 /2d, one 
obtains D^ 1 = AirV^eN 3 / 2 . 

Note that apart from an unimportant numerical pre- 
factor and the logarithm in the denominator of Eq. (74) 
the latter result for t v q coincides with that derived in 
Refs. [37,38,39] for a bulk diffusive metal within the 
framework of a completely different approach, cf., e.g., 
Eq. (81) in [38]. Within that approach local proper- 
ties of the model could not be fully defined. For this 
reason in the corresponding integrals in [37,38,39] we 
could not avoid using an effective high frequency cut- 
off procedure which yields the correct leading depen- 
dence t v o oc D 3 and it only does not allow to recover an 
additional logarithmic dependence on D in (74). Our 
present approach is divergence-free and, hence, it does 
not require any cutoffs. 

We can also add that Eq. (66) also agrees with our 
earlier results [37,38,39] derived for quasi-f d and quasi- 
2d metallic conductors. Provided the transversal size a 
of our array is smaller than d one should set A ~ da for 
2d and A ~ a 2 for Id conductors. Then Eq. (66) yields 
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t v o oc D 2 /\nD and r^o oc D/\nD respectively in 2d 
and Id cases. Up to the factor In D these dependencies 
coincide with ones derived previously, cf., e.g., Eq. (32) 
in [39]. 

Now let us turn to the model (6) of strongly disor- 
dered or granular conductors. In contrast to the sit- 
uation (a), we will assume that the contact between 
dots (grains) is rather poor, and inter-grain electron 
transport may occur only via a limited number of con- 
ducting channels. In this case the average dimension- 
less conductance g can be approximated by some A- 
independent constant g = g c . Substituting g c instead 
of g into Eq. (72) we observe that in the case of strongly 
disordered structures one can expect D oc 1/d. Accord- 
ingly, for Ttpo (66) one finds 

„3 



T V Q = 



327r 2 Ar 2 L> 3 ln(L> c2 /-D)' 



(75) 



where D C 2 again depends on Ec- For Ec ~ e j2d we 
have D~2 — 2-K\/2No/a. Hence, the dependence of t v o 
on D for strongly disordered or granular conductors 
(75) is it qualitatively different from that for suffi- 
ciently clean metals (74). 

One can also roughly estimate the crossover between 
the regimes (a) and (6) by requiring the values of D — 
7«Fd/4 (73) and D = g c /4irNod to be of the same 
order. This condition yields (pFd) 2 ~ 2ivg c /-y, and we 
arrive at the estimate for D at the crossover 



D; 



0.6ft 



(76) 



Here we restored the Planck constant h set equal to 
unity elsewhere in our paper. 

In the next section we will use the above results and 
carry out a detailed comparison between our theory 
and numerous available experimental data for t v q in 
different types of disordered conductors. 



6. Comparison with experiments 

Turning to experiments, it is important to empha- 
size again that low temperature saturation of the elec- 
tron decoherence time has been repeatedly observed 
in numerous experiments and is presently considered 
as firmly established and indisputably existing phe- 
nomenon. At the same time, the physical origin of this 
phenomenon still remains under debate. The key ob- 



servations and remaining controversies are briefly sum- 
marized below. 

(i) The authors [36] have analyzed the values of r^o 
observed in their experiments with open quan- 
tum dots as well in earlier experiments by dif- 
ferent groups [32,33,34,35]. For all 14 samples 
reported in [32,33,34,35,36] the values t v o were 
found to rather closely follow a simple depen- 
dence 



(ii) 



(iv) 



1>0 ~ To- 



(77) 



This approximate scaling was observed within 
the interval of dwell times td of about 3 decades, 
see Fig. 5 in Ref. [36] . To the best of our knowl- 
edge, until now no physical interpretation of this 
observation has been suggested. 
In some of our earlier publications [37,38,45,46] 
we have demonstrated a good quantitative 
agreement between our theoretical predictions 
[37,38] and experimental data for t v o obtained 
for numerous metallic wires and quasi- Id semi- 
conductors. As our theory of dephasing by 
electron-electron interactions [37,38] predicts 
a rather steep increase of r^o with the system 
diffusion coefficient D, e.g. for most metals as 
Ttpo oc D 3 , we can conclude that for a large 
number of disordered conductors r^o strongly 
increases with increasing D. 
In a series of papers, see, e.g., Refs. [26,27,31,47], 
Lin and coworkers analyzed numerous experi- 
mental data for t v q obtained by various groups 
in disordered conductors with D < 10 cm 2 /s 
and observed systematic decrease of t v o with in- 
creasing D. The data could be rather well fitted 
by the dependence t v o oc D~ a with the power 
a > 1. This trend is opposite to one observed in 
less disordered conductors with D > 10 cm 2 /s 
and remained unexplained until now. 
The authors [62] pointed out a disagreement 
between our expressions for t v q [37,38] and the 
data obtained for a number of typically rather 
strongly disordered 2d and 3d structures. In 
some cases this disagreement was argued to be 
as large as 4 to 5 orders of magnitude. In Ref. [66] 
we countered this critique showing, on one hand, 
reasonable agreement for some of the samples in 
question and arguing, on the other hand, that 
our quasiclassical theory [37,38] is applicable 
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merely to weakly disordered conductors. Hence, 
it cannot be used in order to quantitatively de- 
scribe strongly disordered structures like, e.g., 
granular metals, metallic glasses etc. quoted in 
Ref. [62]. Formally eliminating the controversy, 
this our argument, however, did not yet allow to 
clarify the issue of low temperature saturation 
of t v in strongly disordered conductors which 
remained unclear until now. 
(v) Pierre et al. [7] argued that low temperature sat- 
uration of t v in sufficiently clean samples can be 
caused by undetectably small number of mag- 
netic impurities. This idea, however, was not 
supported by the authors [8] who demonstrated 
that at least in their experiments (performed in 
high magnetic fields in order to fully polarize 
any magnetic moments should they exist) the 
observed r^-saturation cannot be due to mag- 
netic impurities. Later the issue was reanalyzed 
in Refs. [12,14,15,16] on the basis of recently de- 
veloped numerical renormalization group (NRG) 
theory of electron scattering by Kondo impuri- 
ties [67]. In samples with implanted magnetic 
impurities a very good agreement between the- 
ory [67] and experiments [14,15,16] was demon- 
strated both above and below the Kondo tem- 
perature Tk- However, at T < O.ITk significant 
deviations from NRG predictions was observed 
and T,f, was found to saturate [14,15,16] in all 
samples both with and without magnetic impu- 
rities. Interpretation of this saturation effect in 
terms of both underscreened and overscreened 
models is problematic [14,16]. 
Let us now analyze the above problems and contro- 
versies point by point within the framework of our the- 
ory of electron-electron interactions. 

We start from the case of quantum dots [32,33,34,35,36] 
observe that our results for t v q (45), (66) scale prac- 
tically linearly with the dwell time td which is essen- 
tially the scaling (77) suggested in Ref. [36]. There is, 
however, one point which requires a comment. 

As it was argued in Sec. 2, one should expect no 
dephasing by electron-electron interactions in single 
quantum dots at any T. Naively one could regard this 
statement as yet one more controversy with experi- 
ments [32,33,34,35,36] where electron dephasing in sin- 
gle dots was clearly observed. At this stage let us recall 
that qualitative arguments in Sec. 2 as well as rigorous 



analysis in Refs. [50,51] remain applicable to single dots 
provided fluctuating voltages drop strictly across the 
two barriers. In realistic quantum dots [32,33,34,35,36] 
fluctuating voltages most likely penetrate inside rather 
than drop only at the edges. Within the framework of 
our model one can easily mimic this situation by in- 
troducing additional scatterers inside the dot in which 
case fluctuating voltages already do dephase (see Sec. 

2). 

For illustration, let us consider a strongly asymmet- 
ric double dot system, i.e. we replace, say, the left bar- 
rier in Fig. 1 by a small quantum dot with the electron 
flight time ta <C td- The left dot will then model a 
barrier of a finite length in a single dot configuration 
of Fig. 1. Applying Eqs. (34), (47) and (48), we evalu- 
ate the WL correction and again arrive at Eqs. (39-40), 
where, however, u = t^Jtrc- Extracting the decoher- 
ence time r^o, we obtain 

Tv " ^ (m/2T R c)V9z - i ( 78 > 

for Tfl > trc and vanishing decoherence rate in the 
opposite limit ta <C trc- Provided the denominator in 
Eq. (78) is not very large, we arrive at Eq. (77). We 
also note that t v o in a single dot (78) is always longer 
than that in a system of two dots (44). 

Alternatively, one can model a dot by a chain of sev- 
eral (N) scatterers with T V (T) defined in Eqs. (66), 
(69), (70). Then at higher temperatures we obtain r v oc 
T~ v with v ranging from 2/3 to 1, as observed in a 
number of dots [32,33,34,35,36]. Substituting qtd by 
fftotT"!?* /N (where g to t are Tjj* are respectively the di- 
mensionless conductance and the dwell time of a com- 
posite dot), bearing in mind that gtot ~ 1 -J- 12 [36] 
and assuming N to be not very large, at low T from 
(66) one finds r^o ~ Td m agreement with experimen- 
tal results presented in Fig. 5 of Ref. [36] . More com- 
plicated configurations of scatterers can also be con- 
sidered with essentially the same results. We conclude 
that our present theory is in a good agreement with 
experimental findings [32,33,34,35,36]. 

Let us now turn to experiments with spatially ex- 
tended conductors. In Fig. 6 we have collected experi- 
mental data for r^o obtained in over 120 metallic sam- 
ples with diffusion coefficients varying by ~ 4 decades, 
from D « 0.3 cm 2 /s to D ~ 350 cm 2 /s. The data were 
taken from about 30 different publications listed in fig- 
ure caption. We see that the the measured values of 
t v o strongly depend on D. Furthermore, this depen- 
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Fig. 6. The low temperature dephasing times observed in 
various experiments: Au [9], [12] (□); Au [8] (A); Au,Ag 
[7] (0); Au Pd [27] (o); [10] (•); Au [15] (V); [24] (B, ten 
points within the box); [31] AU2AI (-4), Sb (<), ScgsAgis 
(►), V3AI (B); [6] (>); [47] CuGeAu (ffl); AuPd-3 and Al 
[70] (♦) Au [25](A) ; . 

dence turns out to be non-monotonous: For relatively 
weakly disordered structures with D > 10 cm 2 /s t v o 
clearly increases with increasing D, while for strongly 
disordered conductors with D < 10 cm 2 /s the oppo- 
site trend takes place. In addition to the data points in 
Fig. 6 we indicate the dependencies t v o(D) (74) and 
(75) for two models (a) and (6) discussed in Sec. 5.2. 

We observe that for D > 10 cm 2 /s the data points 
clearly follow the scaling (74). Practically all data 
points remain within the strip between the two lines 
corresponding to Eq. (74) with 7 = 1 (dashed line) 
and 7 = 0.2 (solid line). On the other hand, for more 
disordered conductors with D < 10 cm 2 /s the data 
are consistent with the scaling (75) obtained within 
the model (6). We would like to emphasize that the- 
oretical curves (74) and (75) are presented in Fig. 6 
without any additional fit parameters except for a ge- 
ometry factor 7 for the first dependence and the value 
g c ss 150 for the second one. This value of g c was esti- 
mated from the crossover condition (76) with D ~ 10 
cm 2 /s and 7 ~ 1. 

Note that quite a few data points with D > 10 
cm 2 /sek correspond to the samples contaminated by 
magnetic impurities. Remarkably, these data points 
also demonstrate - though with somewhat larger scat- 
ter - systematic increase of t v o with increasing D. At 
the same time, for similar values of D the samples with 
higher concentration of magnetic impurities have sys- 



tematically lower T^o than samples with few or no mag- 
netic impurities. These observations indicate that for 
samples with relatively high concentration of magnetic 
impurities both mechanisms of electron-electron in- 
teractions and spin-flip scattering provide substantial 
contributions to t v , being responsible respectively for 
the scaling r^o oc D 3 and for additional non-universal 
shift of the data points downwards. 

In order to carry out more accurate comparison with 
our theory of electron-electron interactions let us now 
leave out the data points taken from the samples with 
high concentration of magnetic impurities. In Fig. 7 we 
selected the data for 37 different metallic wires with 
no or few magnetic impurities and with diffusion co- 
efficients in the range 9 cm 2 /s < D < 300cm 2 /s from 
Refs. [4,5,6,7,8,9,10,11,12,14,15,16,69,70]. Two minor 
adjustments of some data points are in order. Firstly, 
in order to eliminate the uncertainty related to differ- 
ent definitions of t v used by different groups [71], we 
have adjusted 6 data points [6] and 2 data points [10] 
according to the definition of t v used in other works 
[4,5,7,11,12,14,15]. Secondly, 5 data points correspond- 
ing to samples with high diffusion coefficients [7] have 
been adjusted in order to eliminate the temperature 
dependent contribution to t v which remains substan- 
tial for these samples down to the lowest T w 40 mK 
[72]. Both these adjustments can in no way influence 
our conclusions (in fact, non-adjusted data points also 
remain in-between solid and dashed lines in Fig. 7). 

The data in Fig. 7 again clearly support the scal- 
ing of t v o with D (74) due to electron-electron inter- 
actions. Practically all data points are now located in- 
between the solid and dashed lines which indicate the 
dependence (74) with respectively 7 = 0.2 and 7 = 0.5. 
We cannot exclude that the remaining scatter among 
the data points with similar values of D - to a cer- 
tain extent - may be due to relatively small amount 
of magnetic impurities possibly residing in some sam- 
ples. Also minor differences in metal parameters (e.g, 
Fermi velocities) may contribute to this effect. It ap- 
pears, however, that sample-to-sample fluctuations of 
the geometry parameter 7 play the most important 
role. This parameter is defined as a ratio between the 
square root of the inter-dot (inter-grain) contact area 
A 1 ' 2 and the average dot (grain) size V , i.e. 7 ac- 
counts for local properties of the sample which can be 
highly non-universal. Electrons in metallic wires get 
scattered both in the bulk and on the surface. Hence, 
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D (cm7s) 

Fig. 7. The same as in Fig. 6 for 9 cm 2 /s < D < 300cm 2 /s. 
Only the data for the samples with low concentration or 
with no magnetic impurities have been selected: Au-1, Au-3, 
Au-4, Au-6 [4], Au-7 [5], Au-9 [69] and CF-1, CF-2 [8] (■); 
A to F (AuPd) [6] G); Ag(6N)a,b,c,d, Au(6N), Ag(5N)a,b, 
and Cu(6N)a,b [7] (0); Au [9], Aul[12], Agl [16] and Ag2 
[14] (□); 2 samples from [10], still to be included, 2 Au [15] 
(V); D (AuPd), F, H (Ag) [11] (>); AuPd-3 and Al [70](*) 
Au [25] (A). 



e.g., the surface quality and/or details of the wire ge- 
ometry - along with other factors - may significantly 
impact the value of 7. Since r^o (74) depends quite 
strongly on 7, it is by no means surprising that differ- 
ent wires with similar values of D may have decoher- 
ence times differing by few times. 

As an illustration of this point let us consider 4 quasi- 
ld AuPd samples C, D, E, F [6] with nominally iden- 
tical values of D = 15 cm 2 /s. Although these samples 
were fabricated in the same way and measured in one 
experiment [6], their dephasing times were found to dif- 
fer by up to ~ 2 times. This difference can hardly be as- 
cribed to magnetic impurities since the measured level 
of dephasing would require unrealistically high concen- 
tration of such impurities (in the range of 100 ppm) 
and, in addition, rather exotic values of the Kondo 
temperature Tk ■ At the same time, sample-to-sample 
fluctuations of the parameter 7 of only < 25 per cent 
(which is easy to assume given, e.g., somewhat differ- 
ent geometry of the samples) would fully account for 
the above difference of dephasing times. We conclude 
that our theoretical expression (74) is in a good quanti- 
tative agreement with the available experimental data 
for quasi- Id metallic wires with diffusion coefficients in 
the range 9 cm 2 /s < D < 300 cm 2 /s. 



Now let us return to Fig. 6 and consider the data 
for strongly disordered conductors with D < 10 cm 2 /s. 
As we already pointed out, the agreement between the 
data and the dependence (75) predicted within our 
simple model (b) is reasonable, in particular for sam- 
ples with D < 3cm 2 /s. At higher diffusion coefficients 
most of the data points indicate a weaker dependence 
of r^o on D which appears natural in the vicinity of 
the crossover to the dependence (74). The best fit for 
the whole range 0.3 cm /s < D < 10 cm /s is achieved 
with the function t v q oc D~ a with the power a w 
1.5 -j- 2. Further modifications of our - clearly over- 
simplified - model (b) can help to achieve even better 
agreement between theory and experiment. 

Although such modifications can certainly be worked 
out, it is not our aim to do it here. More importantly, 
our analysis of Sec. 5.2 allows to qualitatively under- 
stand and explain seemingly contradicting dependen- 
cies of t v o on D observed in weakly and strongly dis- 
ordered conductors. While the trend "less disorder - 
less decoherence" (74) for sufficiently clean conductors 
is quite obvious, the opposite trend "more disorder - 
less decoherence" in strongly disordered structures re- 
quires a comment. Effectively the latter dependence 
implies that with increasing disorder electrons spend 
more time in the areas with fluctuating in time but 
spatially uniform potentials which do not dephase, as 
we also discussed in Sec. 2. In other words, in this case 
an effective dwell time td in Eq. (66) becomes longer 
with increasing disorder and, hence, the electron deco- 
herence time t v o does so too. 

Is the physical picture of D decreasing with increas- 
ing dot (or grain) size (employed within the model 
(b)) realistic? Although for cleaner conductors the ten- 
dency is usually just the opposite, for strongly disor- 
dered structures increasing resistivity with increasing 
grain size has been observed in a various experiments 
[73,74,75]. In addition, since local conductance fluctu- 
ations increase with increasing disorder, several grains 
can form a cluster with internal inter-grain conduc- 
tances strongly exceeding those at its edges. In this case 
fluctuating potentials remain almost uniform inside the 
whole cluster which will then play a role of an effec- 
tive (bigger) grain/dot. Accordingly, the average vol- 
ume of such "composite dots" V oc 1/8 may grow with 
increasing disorder, electrons will spend more time in 
such bigger dots and, hence, the electron decoherence 
time (66) will increase. 
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The above comparison with experiments con- 
firms that our previous quasiclassical results for t v q 
[37,38,39] are applicable to relatively weakly disor- 
dered structures with D > fO cm 2 /s, while for con- 
ductors with stronger disorder different expressions 
for t v o (e.g., Eq. (75)) should be used. For instance, 
the claimed in Ref. [62] "disagreement" between our 
theory and some experiments by 4 to 5 orders of mag- 
nitude is solely due to misuse of our results [37,38] far 
beyond their applicability range. A glance at Fig. 6 is 
sufficient to realize that any attempt to apply Eq. (74) 
to structures with, e.g., D ~ 0.3 -j- 1 cm 2 /s can eas- 
ily lead to a "disagreement" with the data, say, by 6 
orders of magnitude or so. This observation, however, 
does not imply any real disagreement, rather it indi- 
cates that more general results for t v q, e.g., Eq. (66), 
should be used. The same argument invalidates the 
comparison between our quasiclassical results [37,38] 
and the data for the sample Au-5 [4], see Fig. 8 of Ref. 
[7], or, to a somewhat lesser extent, for the sample C 
of Ref. [11], see Table II of that paper [76]. 

Finally, our analysis allows to rule out scattering 
on magnetic impurities as a cause of low temperature 
saturation of t v . This mechanism can explain neither 
strong and non-trivial dependence of the electron deco- 
herence time on D (see our Figs. 6 and 7 as well as, e.g., 
Fig. 4 in [18], Fig. 1 in [27] and Fig. 5 in [36]) nor even 
the level of dephasing observed in numerous experi- 
ments. E.g., in order to be able to attribute dephasing 
times as short as r^o ^S 10~ 12 s to magnetic impurities 
one needs to assume huge concentration of such impu- 
rities ranging from few hundreds to few thousands ppm 
which appears highly unrealistic, in particular for sys- 
tems like carbon nanotubes, 2DEGs or quantum dots. 
Similar arguments were independently emphasized by 
Lin and coworkers [27,47]. Even in metallic wires with 
high values of D and long dephasing times, at T < 
O.ITk one observes clear saturation of t v [12,14,15] 
which is in a quantitative agreement with our theory of 
electron-electron interactions (see Fig. 7) and is very 
hard to explain otherwise [16]. 

Thus, although electron dephasing due to scattering 
on magnetic impurities is by itself an interesting is- 
sue, its role in low temperature saturation of t v in dis- 
ordered conductors is sometimes strongly overempha- 
sized. Since the latter phenomenon has been repeat- 
edly observed in all types of disordered conductors, the 
physics behind it should most likely be universal and 



fundamental. We believe - and have demonstrated here 
- that it is indeed the case: Zero temperature electron 
decoherence in all types of conductors discussed above 
is caused by electron-electron interactions. 



7. Conclusions 

In this paper we have employed a model of an array 
of quantum dots/scatterers (Fig. 3) which embraces 
various types of disordered conductors and allows to 
study electron transport in the presence of interactions 
within a very general theoretical framework. We have 
non-perturbatively analyzed the impact of electron- 
electron interactions on weak localization for such 
structures with the emphasis put on the interaction- 
induced decoherence of electrons at low temperatures. 
We have formulated a fully self-contained theory free 
of any divergencies and cutoffs which allows to conve- 
niently handle disorder averaging and treat electron 
scattering without residing to quasiclassics. In the 
case of two quantum dots (or three scatterers) it was 
possible to find an exact solution of the problem (Sec. 
4) and to evaluate the WL correction to the system 
conductance practically without approximations. 

With the aid of our approach we have formulated a 
unified description of electron dephasing by Coulomb 
interaction in different structures including (i) weakly 
disordered conductors (e.g., metallic wires with D > 
10 cm 2 /s), (ii) strongly disordered conductors (D < 
10 cm 2 /s) and (iii) metallic quantum dots. We have 
demonstrated that in all these cases at T — > the elec- 
tron decoherence time is determined by the same sim- 
ple formula r^o ~ gi~D/\n(Ec/S). In the case (i) this 
formula yields t v o oc D 3 /lnD and matches with our 
previous quasiclassical results [37,38,39] while in the 
cases (ii) and (iii) it illustrates new physics which was 
not yet explored before. In particular, this formula em- 
phasizes the dependence of t v q on the electron dwell 
time td in single quantum dots [36] and helps to un- 
derstand the (at the first sight counterintuitive) trend 
"more disorder - less decoherence" observed in strongly 
disordered conductors [26,27,31,47]. 

We have carried out a detailed comparison of our 
theoretical predictions with the results of numerous ex- 
periments for the whole scope of structures (i), (ii) and 
(iii) (Sec. 6). In all cases we found a good agreement 
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between theory and experiment which further supports 
our main conclusion that low temperature saturation 
of r v is universally caused by electron-electron inter- 
actions. 
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